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Abstract 

In central collisions at relativistic heavy ion colliders like the Relativistic 
Heavy Ion Collider RHIC/Brookhaven and the Large Hadron Collider LHC 
(in its heavy ion mode) at CERN/Geneva, one aims at detecting a new form of 
hadronic matter — the Quark Gluon Plasma. It is the purpose of this review 
to discuss a complementary aspect of these collisions, the very peripheral ones. 
Due to coherence, there are strong electromagnetic fields of short duration in 
such collisions. They give rise to photon-photon and photon-nucleus collisions 
with high flux up to an invariant mass region hitherto unexplored experimen- 
tally. After a general survey photon-photon luminosities in relativistic heavy 
ion collisions are discussed. Special care is taken to include the effects of 
strong interactions and nuclear size. Then photon-photon physics at various 
77-invariant mass scales is discussed. The region of several GeV, relevant for 
RHIC is dominated by QCD phenomena (meson and vector meson pair pro- 
duction). Invariant masses of up to about 100 GeV can be reached at LHC, 
and the potential for new physics is discussed. Photonuclear reactions and 
other important background effects, especially diffractive processes are also 
discussed. A special chapter is devoted to lepton-pair production, especially 
electron-positron pair production; due to the strong fields new phenomena, 
especially multiple e + e~ pair production, will occur there. 
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Q 2 = ^ + Ql (2) 



1 Introduction 

The parton picture is very useful to study scattering processes at very high energies. 
In this model the scattering is described as an incoherent superposition of the 
scattering of the various constituents. For example, nuclei consist of nucleons which 
in turn consist of quarks and gluons, photons consist of lepton pairs, electrons consist 
of photons, etc. It is the subject of this topical review to discuss that relativistic 
nuclei have photons as an important constituent, especially for low enough virtuality 
Q 2 = — q 2 > of the photon. This is due to the coherent action of all the charges 
in the nucleus. The virtuality of the photon is related to the size R of the nucleus 
by 

Q 2 <i/i? 2 , (l) 

the condition for coherence. The radius of a nucleus is given approximately by 
R = 1.2 fm A 1 / 3 , where A is the nucleoli number. From the kinematics of the 
process one has 

u? 
7 2 

Due to the coherence condition the maximum energy of the quasireal photon is 
therefore given by 

bJmax ~ ~fZi (3) 

H 

and the maximum value of the perpendicular component is given by 

Q± Z (4) 

We define the ratio x = uj/E, where E denotes the energy of the nucleus E = M^jA 
and Mn is the nucleon mass. It is therefore smaller than 

max ~ RM N A ~ R ' [ ' 

where Ac (A) is the Compton wave length of the ion. Here and also throughout the 
rest of the paper we use natural units, setting h = c = 1. 

The collisions of e + and e~ has been the traditional way to study 77-collisions. 
Similarly photon- photon collisions can also be observed in hadron-hadron collisions. 
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Since the photon number scales with Z" 1 (Z being the charge number of the nucleus) 
such effects can be particularly large. Of course, the strong interaction of the two 
nuclei has to be taken into consideration. 



z 




x=co/E 



Figure 1: A fast moving nucleus with charge Ze is surrounded by a strong elec- 
tromagnetic field. This can be viewed as a cloud of virtual photons. These pho- 
tons can often be considered as real. They are called equivalent or quasireal pho- 
tons. The ratio of the photon energy uj and the incident beam energy E is de- 
noted by x = uj / E . Its maximal value is restricted by the coherence condition to 
x < X C (A)/R re 0.175/A 4 / 3 , that is, x < 10 3 for Ca ions and x < 10 4 for Pb 
ions. 




Figure 2: Two fast moving electrically charged objects are an abundant source of 
(quasireal) photons. They can collide with each other and with the other nucleus. 
For peripheral collisions with impact parameters b > 2R, this is useful for photon- 
photon as well as photon-nucleus collisions. 

The equivalent photon flux present in medium and high energy nuclear collisions 
is very high, and has found many useful applications in nuclear physics [Q , nuclear 
astrophysics particle physics Q (sometimes called the "Primakoff effect"), as 

well as, atomic physics ||. It is the main purpose of this review to discuss the physics 
of photon-photon and photon-hadron (nucleus) collisions in high energy heavy ion 
collisions. With the construction of the "Relativistic Heavy Ion Collider" (RHIC) 
and the "Large Hadron Collider" (LHC) scheduled for 1999 and for 2004/2008, 
respectively, one will be able to investigate such collisions experimentally. The 
main purpose of these heavy ion colliders is the formation and detection of the 
quark-gluon-plasma, a new form of highly excited dense hadronic matter. Such 
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a state of matter will be created in central collisions. The present interest is in 
the "very peripheral (distant) collisions" , where the nuclei do not interact strongly 
with each other. From this point of view, grazing collisions and central collisions 
are considered as a background. It is needless to say that this "background" can 
also be interesting physics of its own. 

The equivalent photon spectrum extends up to several GeV at RHIC energies 
(7 « 100) and up to about 100 GeV at LHC energies (7 « 3000), see Eq. (g). 
Therefore the range of invariant masses M 77 at RHIC will be up to about the mass 
of the rjc, at LHC it will extend into an invariant mass range hitherto unexplored. 

We discuss the equivalent photon spectra of strongly interacting particles, from 
which 77- luminosities are obtained. Due to the coherence effect, the corresponding 
77-luminosity is very high. In addition higher order and inelastic processes, which 
may occur in heavy ion collisions are discussed. Then the various possibilities for 
77-physics in the different invariant mass regions will be explored. A special case of 
77-physics is lepton pair production, in particular the e + e~ pair production. Since 
the equivalent photon approximation fails in some regions of phase space, this case 
is discussed separately. Due to its strong-field aspects, the production of multiple 
e + e~-pairs is of basic interest. Electron-positron pair creation is also of practical 
interest due to the possibility of capturing an electron in a K,L,. . . shell. This 
changes the charge state of the ion and leads to a beam loss and thus to a decrease 
in the luminosity. The very large cross section for the production of e + e~ pairs 
(sometimes called QED-electrons) is also an important background for detectors. 

Relativistic heavy ion collisions have been suggested as a general tool for two 
photon physics about a decade ago. Yet the study of a special case, the produc- 
tion of e + e _ pairs in nucleus-nucleus collisions, goes back to the work of Landau 
and Lifschitz in 1934 || (In those days, of course, one thought more about high 
energy cosmic ray nuclei than relativistic heavy ion colliders). In the meantime 
the importance of this process has become very clear, and many studies followed, 
e.g., This subject will be dealt with in detail in Sec. |B| and ^, where also 

recent experimental results will be mentioned. 

The general possibilities and characteristic features of two-photon physics in rel- 
ativistic heavy ion collisions have been discussed injBj. The possibility to produce 
a Higgs boson via 77-fusion was suggested in (h], [yj. In these papers the effect 
of strong absorption in heavy ion collisions was not taken into account. This ab- 
sorption is a feature, which is quite different from the two-photon physics at e + e~ 
colliders. The problem of taking strong interactions into account was solved by using 
impact parameter space methods in [^-[Tif . Thus the calculation of 77-luminosities 
in heavy ion collisions is put on a firm basis and rather definite conclusions were 
reached by many groups working in the field [fl5|| ; for a recent review containing 
further references see |fL6fl . Subsequent studies — to be described in detail in this 
review — revealed in a clear way that the theoretical situation is basically under- 
stood. This opens the way for many interesting applications. Up to now hadron- 
hadron collisions have not been used for two-photon physics. An exception can be 
found in |jl7). There the production of ^+ [iT pairs at the ISR was observed. The 
special class of events was selected, where no hadrons are seen associated with the 
muon pair in a large solid angle vertex detector. In this way one makes sure that the 
hadrons do not interact strongly with each other, i.e., one is dealing with peripheral 
collisions (with impact parameters b > 2R); the photon-photon collisions manifest 
themselves as "silent events" . 

We feel that this is a very good basis for planning concrete experiments, as it 
is done at RHIC ]Tg|-|2l|] and LHC p^-p4|. This review aims at giving the main 
physical ideas and provide the key formula and results. Details can be found in the 
literature. A few new results will also be presented. But the main emphasis is to 
discuss the principle ideas and results in the field. We hope that this review will 
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further stimulate future investigations. It is appropriate to recall that RHIC will 
start operating in 1999, only a year from now. 



2 General survey of peripheral collisions 

Let us first discuss the importance of electromagnetic interactions in peripheral 
collisions for the case of elastic scattering. The strength of the Coulomb interaction 
is measured by the Coulomb (or Sommerfeld) parameter rj which is given in terms 
of the nuclear charges Z\ and Zi by 

ZxZ 2 e 2 

??=— ? ■ (6) 

nv 

For (ultra)relativistic collisions we have v ss c and thus rj « Z\Z-2,jYil . There- 
fore for pp-collisions we have always rj <C 1 and the Born approximation, that is, 
one-photon exchange, is applicable. For this case, i.e., for 77 w 1/137, elastic scat- 
tering is reviewed, e.g., in p5| . (Experimental result for pp scattering at y/s = 546 
and 1800 GeV at Fermilab are given in (2^ |2S|]). There is a superposition of nu- 
clear and Coulomb amplitudes and the elastic differential scattering cross section 
is divided into three distinct regions, separated by the value of the square of the 
momentum transfer t. It is defined as usual as t — (pi — Pf) 2 , which is negative 
(space-like momentum transfer) in our metric used. For \t\ <C \t\i n t Coulomb scat- 
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tcring dominates, for \t\ 3> \t\i n u nuclear scattering dominates and for t\ ~ \t 
there is Coulomb-nuclear interference. \t\i nt is given by (see Eq. (3.13) of |2q|, where 
the factor of Z\Zi is added for the scattering of two nuclei with charge Z\ ,Z% instead 
of two protons): 

|*|*- « (7) 
The nuclear elastic scattering amplitude is usually parameterized as 

F n = iE±p!l e Bt/^ (8) 

where a to t is the total cross section, p = and B is a slope parameter 

related to the size of the hadron. In this normalization the differential cross section 
is given by 

^ = \F n f. (9) 

On the other hand for (very) heavy ions we have 77 3> 1 and semi-classical 
methods are appropriate to deal with elastic scattering (see also |p9| ). For 77 ^> 1 the 
Coulomb interaction is very strong and the Born approximation is no longer valid. 
(In an analogy with optics there is now Fraunhofer diffraction instead of Fresnel 
diffraction). Instead one should use a Glauber approximation, where the Coulomb 
interaction is taken into account to all orders. One can say that many photons are 
typically exchanged in the clastic collisions (in contrast to the Born approximation 
relevant for the pp case, where one photon exchange is sufficient). One can now 
integrate over the impact parameter b using the saddle point approximation. Thus 
one recovers the classical picture of scattering. The particles move essentially on a 
straight line with a constant velocity and an impact parameter b. For Z\ = Z 2 = Z 
at grazing impact parameter b — 2R the momentum transfer is given by (see, 
e.g.,©) 

/A ,2 {2Z 2 e 2 \ 2 Z 4 e 4 
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This momentum transfer determines the scattering angle 6 — Ap/p, where p is the 
momentum of the particle. 

Let us compare this quantity with the corresponding one for pp-scattering, 
Eq. (0). With a tot ~ 8?ri? 2 one has 

\t\int = (11) 

to be compared to 

n/Z 2 

\t\com = -j^{aZ 2 ). (12) 

It seems interesting to note that |i|cotii is Z 2 a times the corresponding quantity 
\t\i n t. Evidently the corresponding scattering angles are exceedingly small. For 
i] w l,i.e., for Z w 12, the quantities \t\i nt and |i|co«z (Eq. jll| ) and (|l^)) are about 
equal, as should be the case. For such values of Z the change from Fresnel to 
Fraunhofer scattering takes place. A more detailed discussion is given in Ch. 5.3.4 
of flj, based on the work of §§. 

The pp cross section is rising with energy (see, e.g., p. 193 of p3[]). At LHC 
energies it will be of the order of 80 mb. What does this mean for p-A and A- 
A cross section? A discussion is given in p^ ]. Experimental information can be 
obtained from cosmic ray data for p-A reactions. Calculation for cross section for 
A-A data often make use of a density-folding approach. This approach can be 
justified starting from Glauber theory fl35| , |36| ■ One uses the thickness function 
Ta(6) of each nucleus, which is defined as the projection of the density along the 
beam axis: 

dz n A Ub 2 + z 2 ) . (13) 

-oo 

Here the nuclear density n^(r) is normalized to J d 3 r riA{r) = A. From this one 
gets the "overlap function" given by 

Tab (b) = J d 2 b x J d 2 b 2 t(b- h - b 2 )T A (b 1 )T B (b 2 ), (14) 

where & is the impact parameter between the two ions (see also Fig. ||). 

t(s) describes the finite range of the nucleon-nucleon interaction. It is propor- 
tional to the profile function T(b) in the Glauber theory, but it is normalized to 
J d 2 s t(s) = 1. Therefore we have J d 2 bTAB{b) = AB. For high energies the NN 
cross section rises beyond the geometrical size of the protons. The finite range of 
the nucleon-nucleon interaction is important then. A similar situation occurs also in 
nuclear physics at lower energies, see, e.g., Figs. 2 and 3 of |37|], where the influence 
of the finite range of the NN interaction is clearly seen on the total cross-section for 
scattering of Li and Be isotopes. 

If the nucleon-nucleon cross section is almost purely absorptive, the scattering 
amplitude is almost imaginary and one can get t(b) from the elastic differential cross 
section 

fd 2 P± exp(ip ± s) x /5|^ 
t(s) = ^_ V \ (15) 



(2n) 2 



da 



d- l 



p±=o 



with gg^- the elastic nucleon-nucleon cross section. Often t(b) is approximated by 

an exponential function of the form t(b) ~ exp(— (6/&o) 2 ), where bo is the range of 
the interaction. 
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If one can neglect the finite range of the interaction compared to the geometrical 
size of the nuclei, the "thickness function" can be simplified: 

Tab (6) = / d 2 &i J d 2 b 2 5 {2 \b - 61 - b 2 )T A (b 1 )T B {b2) (16) 

= J d 2 hT A (bx)T B (b-bi). (17) 

The inelastic scattering cross section of the nucleus- nucleus collision is given by 

<TAB,inel = / ^b I 1 _ ex P (~ & N nTab (b))] , (18) 



where cat at is the total nucleon-nucleon cross section. This integration over b allows 
the integrand to be interpreted as a probability for the two ions to interact. This 
probability is about 1 for small impact parameters, where the nuclei overlap, it 
drops to zero for large impact parameters. The width of the area, where it falls 
off, is given by the surface diffuseness of the nuclear densities and the range of the 
nucleon-nucleon interaction as given by t. Similarly the elastic cross section is given 
by 

v A B,el = J d 2 b [1 - exp (-a NN T AB (b)/2)} 2 . (19) 

Using results from Regge theory (see, e.g., ||||3£j) one can estimate the range 
bo of the interaction by 

b Q = 2yjB , P P + \n (s/s )a' p . (20) 



Using for B , P p ~ 2.4GcV~ 2 and for a' P « 0.25GeV~ 2 , as found in @,||], one gets 
for LHC energies (s = (2 x 7TeV) 2 ,so = lGeV 2 ) bo ~ 1 fm, which is small compared 
to the nuclear radius. 

As we will see in the next section, the photon-photon luminosity depends also on 
the probability for the two nuclei to interact with each other. Only those processes, 
where the ions do not interact, are useful for photon-photon physics. From our 
discussion here, we can conclude that the effect of the increasing range of the NN 
interaction will only be of some importance at the very high end of the invariant 
77-masses for a quantitative determination of the photon-photon luminosity. 



a) side view b) head on view 




Figure 3: The probability of the nucleus to interact in central collisions can be 
described within the p-p folding approach to the Glauber theory. The parameters b, 
b\ and bi used in the text are explained. 
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3 From impact-parameter dependent equivalent 
photon spectra to 77-luminosities 

Photon-photon collisions have been studied extensively at e + e~ colliders. The the- 
oretical framework is reviewed, e.g., in EG]. The basic graph for the two-photon 
process in ion-ion collisions is shown in Fig. 14. Two virtual (space-like) photons col- 
lide to form a final state /. In the equivalent photon approximation it is assumed 
that the square of the 4- momentum of the virtual photons is small, i.e., q\ ~ q 2 « 
and the photons can be treated as quasircal. In this case the 77-production is 
factorized into an elementary cross section for the process 7 + 7 — > / (with real 
photons, i.e., q 2 — 0) and a 77-luminosity function. In contrast to the pointlikc 
elementary electrons (positrons), nuclei are extended, strongly interacting objects 
with internal structure. This gives rise to modifications in the theoretical treatment 
of two photon processes. 

3.1 Elastic vertices 

The emission of a photon depends on the (elastic) form factor. Often a Gaussian 
form factor or one of a homogeneous charged sphere is used. The typical behavior 
of a form factor is 

f{q } ~\ for| 9 2 |»^ ■ ^ 

For low I q 2 1 all the protons inside the nucleus act coherently, whereas for | q 2 1 ^> 1 /R 2 
the form factor is very small, close to 0. For a medium size nucleus with, say, R = 5 
fm, the limiting Q 2 = -q 2 = l/R 2 is given by Q 2 = (40MeV) 2 = 1.6 x 10~ 3 GeV 2 . 
Apart from e + e~ (and to a certain extent also /i + /i~) pair production — which will 
be discussed separately below — this scale is much smaller than typical scales in 
the two-photon processes. Therefore the virtual photons in relativistic heavy ion 
collisions can be treated as quasireal. This is a limitation as compared to e + e~ 
collisions, where the two-photon processes can also be studied as a function of the 
corresponding masses q\ and q 2 of the exchanged photon ( "tagged mode" ) . 




Figure 4: The general Feynman diagram of photon-photon processes in heavy ion 
collisions: Two (virtual) photons fuse in a charged particle collision into a final 
system f. 

As was discussed already in the previous section, relativistic heavy ions interact 
strongly when the impact parameter is smaller than the sum of the radii of the 
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Figure 5: "Elastic photon emission": In order for the photon to interact coherently 
with the whole nucleus, its virtuality Q 2 is restricted to Q 2 < 1/R 2 . In the calcu- 
lation this is incorporated by the elastic form factor f(q). This is important for a 
treatment, where plane waves are used. In the semiclassical (or Glauber) method, 
the detailed form of f(q 2 ) is not relevant for collisions with b> R (see text). 




Figure 6: "Inelastic photon emission": As a nucleus is a rather weakly bound system, 
photon emission can lead to its breakup or excitation, especially for Q 2 > 1/R 2 . Also 
for this incoherent case an equivalent photon number, similar to the coherent case, 
can be defined. 
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two nuclei. In such cases 77-processes are still present and are a background that 
has to be considered in central collisions. In order to study "clean" photon-photon 
events however, they have to be eliminated in the calculation of photon-photon 
luminosities as the particle production due to the strong interaction dominates. In 
the usual treatment of photon-photon processes in e + e~ collisions plane waves are 
used and there is no direct information on the impact parameter. For heavy ion 
collisions on the other hand it is very appropriate to introduce impact parameter 
dependent equivalent photon numbers. They have been widely discussed in the 
literature (see e.g. §,|o[|l|) 

The equivalent photon spectrum corresponding to a point charge Ze, moving 
with a velocity v at impact parameter b is given by 



N(u,b) 



Z 2 a 1 



(0 



Kl(x) + ^K 2 (x) 



(22) 



where K n {x) are the modified Bessel Functions (MacDonald Functions) jl2) and 
x = Then one obtains the probability for a certain electromagnetic process to 
occur in terms of the same process generated by an equivalent pulse of light as 



P(b) = J ^JV(o>,&)(7>). 



(23) 



Possible modifications of N(co, b) due to an extended spherically symmetric charge 
distribution arc given in |Q (see also Eq. ( |l9|) below). It should be noted that 
Eq. ( f22"l) also describes the equivalent photon spectrum of an extended charge dis- 
tribution, such as a nucleus, as long as b is larger than the extension of the object. 
This is due to the fact that the electric field of a spherically symmetric system 
depends only on the total charge, which is inside it. As one often wants to avoid 
also final state interaction between the produced system and the nuclei, one has to 
restrict oneself to bi > Ri and therefore the form factor is not very important. For 
inelasti c ve rtices a photon number N(aj,b) can also be defined, as will be discussed 
in Sec. 3.2 below. 




Figure 7: View of the collision perpendicular to the beam direction: The electric 
field vector points along the direction of the individual impact parameter. 

As the term x 2 [Kf(x) + 1 / 'f 2 Kq{x)\ in Eq. ( p^ ) can be roughly approximated 
as 1 for x < 1 and for x > 1, so that the equivalent photon number N(ui,b) 
is almost a constant up to a maximum u> m ax — l/b (x = 1). By integrating the 
photon spectrum (Eq. (p2|)) over b from a minimum value of R m in up to infinity 
(where essentially only impact parameter up to b max j/uj contribute, compare 
with Eq. (pi)), one can define an equivalent photon number n{uo). This integral can 
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be carried out analytically and is given by p],pO[ 

2 „o /c\2 " 



n(w) = / d 2 bN(Lu,b) = -Z\a (^j 



(24) 



^2 
2c 2 

where the argument of the modified Bessel functions is £ = " Rmin . The cross section 
for a certain electromagnetic process is then 

/ddJ 
— n(u)a^(ui). (25) 
UJ 

Using the approximation above for the MacDonald functions, we get an approxi- 
mated form, which is quite reasonable and is useful for estimates: 

( \ ~ 2Z2a i 1 (oc\ 
n{uj) fa In — . (26) 

7T OJHmin 

The photon-photon production cross-section is obtained in a similar factorized 
form, by folding the corresponding equivalent photon spectra of the two colliding 
heavy ions 13, L4| (for polarization effects see they are neglected here) 



^1 f ^Ip^^y = V4^), (27) 



with 



1-2 71 



poo i>oo />z 

F{u x ,U2) = 2tt / hdh / M&2 / 
JRi Jr 2 Jo 
xJV(wi, 6i)JV(w2, 6 2 )6 (6 2 + 6* _ 2M2 cos0 _ i? 2 uto// ) , (28) 

(Rcutoff = -Ri + iia)- This can also be rewritten in terms of the invariant mass 
IU 77 = V4wiw 2 and the rapidity Y = l/21n((P + P z )/(Po - P*)) = 1/2 ln(wi/cj 2 ) 
as: 

^rr^S^y-rrCWn), (29) 
with 

dW^dY - \ 2 6 ' 2 6 J' W 

Here energy and momentum in the beam direction are denoted by Pq and P z . The 
transverse momentum is of the order of P± < 1/R and is neglected here. The 
transverse momentum distribution is calculated in EJ . 

In j44| and |45|] this intuitively plausible formula is derived ab initio, starting 
from the assumption that the two ions move on a straight line with impact parameter 
b. Eqs. ( p7| ) and ( |2^ ) are the basic formulae for 77-physics in relativistic heavy- 
ion collisions. The advantage of heavy nuclei is seen in the coherence factor Z\Z^ 
contained in Eqs. (||)-©. 

Let us make a few remarks: In Eq. a sharp cut-off in impact parameter 
space at b m m = Rcutoff is introduced. There is some ambiguity in the numerical 
value of Rcutoff'- As was discussed in Sec. || the total cross section for nucleon- 
nucleon collisions is rising with energy. To a certain extent this will also affect the 
values of R cu toff- They are not just the (energy independent) radii of the nuclei, 
but they rather describe the probability of nuclear interactions between the heavy 
ions. A more realistic calculation can be done by replacing this sharp cutoff with 
a smooth one, using the overlap function of Eq. ( [H] ) and (|l4|). The O-function in 
Eq. (p|h has to be replaced by 

P(b) = 1 - cxp (-<J NN T AB (b)) , (31) 
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with b = b\ — 62 ■ Comparing this refined model with the one with a sharp cutoff, 
significant deviations are only present at the very upper end of the invariant mass 
range. Only the smallest impact parameter contribute significantly to these events, 
therefore a sensitivity to the cutoff is expected. All other (smaller) invariant masses 
are not very sensitive to the details of this cutoff. 
Integrating over all Y the cross section is 



a c = J 



dL 

dwT 



(32) 



77 



For symmetric collisions (R± = R% = R) we can write the luminosity in terms of a 
universal function /(r) as 



dL/dW, 



77 



Z 4 a 2 1 

TT 4 W2. 



■fir), 



(33) 



with t — Wj-yR/j [jl4|. For large values of 7, f(r) is given by 

/>00 /-OO />2"7T 

/(t) = 7TT 2 / U\dU\ I U2dU2 I d(f) O {u\ + u\ — 2U\U2 COS (j) — 4) 

Jl Jl Jo 

"+°° ^ ^ 

;u 2 e ). 



/-|-00 
dYK\C-u,e Y )KlC^ 



(34) 



The function /(t) is shown in Fig. g, a useful parameterization of it was studied 
in [HI. 




Figure 8: Plot of the universal function fir) (see Eq. 
W^R/j. 



As a function of Y, the luminosity d 2 L j dW^dY has a Gaussian shape with the 
maximum at Y = 0. The width is approximately given by AY = 2 In [(27) / (R W 71 )) . 
Depending on the experimental situation additional cuts in the allowed Y range are 
needed. 
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Ion 


Z,A 




7 


L AA (cm 2 s L ) 


LHC 


Pb-Pb 


82,208 


1148 TeV 


2950 


1 x 10 2e 


Ca-Ca 


20,40 


280 TeV 


3750 


4 x 10 30 


P-P 


1,1 


14 TeV 


7450 


10 29 - 10 31 


RHIC 


Au-Au 


79,197 


20 TeV 


106 


2 x 10 2b 


I-I 


53,127 


13 TeV 


111 


3 x 10 27 


Cu-Cu 


29,63 


7.2 TeV 


122 


1 x 10 28 


Si-Si 


14,28 


3.5 TeV 


133 


4 x 10 28 


0-0 


8,16 


2.0 TeV 


133 


1 x 10 29 


p-p 


1,1 


250 GeV 


266 


1 x 10 31 



Table 1 : Parameters for different ion species at RHIC and LHC. In the entries we 
give the total invariant mass y/s of the system, the Lorentz factor 7 and the beam 
luminosity. Parameters are taken from 



3.2 Inelastic Vertices and Higher Order Corrections 

Heavy ions are complex objects unlike the pointlike, structureless electrons. The 
effects due to the finite size (and the nuclear interactions between them) has been 
considered by using an impact parameter approach. The additional processes com- 
ing from elastic nuclear interactions (diffractive processes, Pomeron interactions) 
will be considered separately in Sec. ^ and [| Here we want to discuss mainly 
additional electromagnetic processes that can occur in distant collisions. 

Especially for very heavy ions — like Pb — Za is no longer small. Therefore 
processes with more than one photon-exchange are important. In addition to the 
77-process one can have electromagnetic dissociation of the ions. Furthermore the 
ions can also be excited due to the emission of the photon. This excitation can lead 
to an excited nucleus, or to the breakup of the nucleus, when the proton is knocked 
out of the nucleus due to the photon emission. At even higher Q 2 the photon can 
be emitted also from the quarks contained in the protons. 

In Sec. |J we will see that the electromagnetic dissociation of the ions is an 
important process. This cross section is often so large that it can also occur in 
addition to the 77-process. 

In the impact parameter dependent approach the probability of several pro- 
cesses to occur in one collision can be calculated by assuming that the processes are 
independent of each other; their probabilities have to multiplied: 

P fA . (b) = P^ f (b) P lA ^ A , (b). (35) 

Therefore the processes given in Fig. ^ are effectively included in this approach. 
This is a good approximation as long as the nucleus is not disturbed substantially. 
Integrating from b = 2R up to infinity, one obtains the cross section for 77 — > / 
fusion accompanied by a specific jA — > A* interaction 

/>oo 

a fA * = 2tt / bdb P^ } {b) P 7X _» A .(6). (36) 

J2R 

Since J2 A * P-yA->A* = 1 (where we sum over all possible states A* of the nucleus 
including the ground state), we have 

^ct m . = CT/) (37) 

A* 
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Figure 9: Due to the large charge of the ions, electromagnetic excitation in addition 
to the photon-photon process may occur. 



that is, the same as the cross section of the 77 — ► / process alone, without any 
higher order processes, as given above. 

One can make an estimate of the importance of the higher order processes in 
the following way: The luminosity especially at high invariant masses peaks at 
b « b m in = 2R and one can therefore make the approximation 

a ftA m P A {2R)a 71 ^ f , (38) 

for the cross section <t/,a, that is 77 — » / with the nucleus in the ground state in 
the final state and 

a ftA ^(l-P A (2R))a^ f (39) 

for the nucleus to be excited. The impact parameter dependent probability for the 
excitation of the nucleus can be found in p3||4§|] . They are strongly dependent 
on A. Especially important (see also Sec. 0) is the giant dipole resonance (GDR), 
a highly coherent excited state, which can be interpreted as the movement of all 
neutrons against all protons. The probability of such an excitation of a nucleus 
with charge Z2 and nucleon number A2 (N2 = A2 — Z2) by a nucleus with charge 
number Z\ is given in |jj as 

Pgdr(P) = 1 - cxp(-S/6 2 ), (40) 

with 

S = 5.45xl0- 5 ^|^/m 2 , (41) 

where the Thomas- Reiche-Kuhn (TRK) sum rule was used and the position of the 
GDR is given by SOA' 1 / 3 MeV. For Pb-Pb we get a value of S = (10.4/m) 2 . A 
parameterization which includes also all additional contributions (quasideuteron 



absorption, nucleon excitation, etc.) was given in 49 1; they get a value of S = 
(17.4/m) 2 . For the system Ca-Ca we get S = (0.86/m) 2 using Eq. (glj). 

From this estimate one finds that about 75% of the photon-photon events in 
Pb-Pb collisions are accompanied by an excitation, about 40% leading to the GDR. 
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For Ca-Ca these are less than 2% of all events. As each of the ions can be excited, 
we get a total probability for excitation of at least one of them of about 2% for 
Ca-Ca and of about 95% for Pb-Pb (65% from GDR). These effects are therefore 
dominant for Pb-Pb, but almost unimportant for Ca-Ca. 

A detailed calculation of the cross section using the full b dependence is shown 
in Fig. |l(]. It shows that the excitation probability is lower than the estimate given 
above. 




j -5 I 1 1 1 1 

50 100 150 200 

W(GeV) 



Figure 10: For Pb-Pb collisions a substantial portion of the total luminosity (o~f, 
solid curve) comes from events with are accompanied by the excitation of at least 
one ion (a a*, double-dashed curve) The dominant contribution is to the giant dipole 
resonance (oqdr, dotted curve). For Ca-Ca collisions this is no longer very im- 
portant. Shown are full calculation for the LHC (7 ~ 2950 for Pb-Pb, 3750 for 
Ca-Ca). 

The GDR excitation is followed most of the time by neutron evaporation. Also 
other photon induced reactions predominantly lead to the emission of individual 
nucleons. This emission of relatively low energy nucleons in the nucleus rest frame 
leads to high energetic neutrons (protons) with energies of about 3 TeV (LHC) or 100 
GeV (RHIC). The neutrons can possibly be detected in a zero degree calorimeter. 

In f30f it was proposed to use the mutual emission of neutrons from both nuclei as 
a measure of the beam luminosity at RHIC. Using the coincidence of two neutrons in 
the very forward and backward direction other sources of neutrons can be suppressed 
effectively Since the photonuclear processes are large and well understood they can 
lead to a good determination of the luminosity. The authors of ]5(| estimate to 
be able to determine the luminosity to about 5%. It seems interesting to note 
that the A dependence of the excitation cross section with two photons is given 
approximately by 10 _9 A 6 mb Jl3], the one-photon exchange cross section is given 
approximately by 10~ 5 A 2 mb ]51| , see Fig. [ll]. So for nuclei heavier than about C, 
the two-photon mechanism is dominant over the one-photon mechanism. 

In the calculation of the luminosity we were always assuming that both nuclei 
remain in their ground state. But nuclei are weakly bound composite systems and 
it is possible that the photon emission leads to their excitation (see also Fig. |i~2| ). 
We distinguish two different types: those leading to an excited nucleus with a 
well defined excitation and the incoherent photon emission from individual protons 
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GDR 



GDR 



Figure 11: The mutual excitation process (a), where both ions are excited due to 
one-photon exchange, becomes less important at larger A compared to the second 
order process (b). 



(or even quarks within the protons), which are best treated as inclusive processes, 
summing over all excitation energies. 

An equation for the inelastic photon-emission to a discrete state, was derived 
in p2| (see also Q]) using plane waves and therefore not subtracting the central 
collisions. It was applied to nuclear excitation, as well as to the case of the proton- 
Delta transition p3]. The equivalent photon number can be expressed in terms of 
the structure functions C and D of the general hadron tensor 



MiM, 



q^q v ~ 


c + 






- pv _ ^ 


1 




1 




1 



D, (42) 



where denotes the nuclear four-current. One obtains 



n(w) 



(-2C + q ± /w 2 P 2 D)m 2 2 
(27r) 3 2 J BP(q 2 ) 2 ' 



(43) 



where E and P are energy and momentum of the nucleus. Whereas in the elastic 
case q 2 was given by —(uj 2 /"/ 2 + q\), it is here s» — — [~ + 2A^j , where A is the 

excitation energy. Using the Goldhaber- Teller model, the transition to the GDR 
was found to be very small, below 1% of the elastic contribution. Using the elastic 
structure functions of the proton in its usual dipole form (see, e.g., p0| ) the elastic 
proton equivalent photon number was obtained analytically in p4| . Quite similarly 
the equivalent photon number corresponding to the p — A transition was obtained 
analytically in p3| using the structure functions of [d5[ as 



^p^a(w) 



with 



*9 / J; , \ 2 

a /i / to + m 



47r 9to 2 \ 2to 
3t 2 



In 
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2t, 



A 2 + / . /' A 2 1 / . 

A 8 



2(A 2 + 3(A 2 + t m i n j J 3(A 2 + tmin) 

uJ 1 (m* — m) c<j(m* — m) 



7 2 



(44) 



(45) 



7 7 

and where to* = 1232 MeV is the mass of the A, m = 938 MeV the proton mass, 
A 2 = 0.71 GeV 2 and \i* — 9.42. For not too large photon energy u> it is given by a 
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constant 

f ,\~ a V* 2 ( m*+m \ 4 A 2 

n ^ )K ^^\r^r) t- (46) 

This is an effect of the order of 10% SI. 




Figure 12: /n photon-photon collisions, one (Fig. (a)) and even both (Fig. (b)) 
nuclear vertices can be inelastic, leading to excited nuclei. 

As we have seen above, the knowledge of impact parameter dependent equiva- 
lent photon numbers is very helpful, as it allows in a direct way to take the strong 
absorption into account (Eq. (|2^)). The equivalent photon spectrum corresponding 
to a fast moving point particle is given by Eq. J22|) . This result can now be general- 
ized to arbitrary charge-current distributions (also quantum-mechanical transition 
currents in the framework of the Glauber theory). 

The equivalent photon spectrum corresponding to a spherically symmetric 
charge distribution p(r) moving with velocity v at an impact parameter b is derived 
in [J43|| - The Fourier transformation of this charge distribution, with J d 3 rp(r) = Z, 
is given by 

f(k 2 ) = ZF{k 2 ) = { d 3 rexp{-ikf)p(r). (47) 



For 7> 1 we only need the perpendicular component of the electric and mag- 
netic fields, E± and B±. One finds E| 

f d 2 k i 7,p - -> - 

^ = / 77T^-^f(k 2 )exp(-tk ± b)k ± , (48) 
J (27r) z vk z 

where k 2 w — k\ — (ui/j) 2 . It can easily be seen, that E± || b and one obtains the 
equivalent photon number as 



N ^ = — U v 2 



, ft x'+u' \ 

duu 2 J 1 (u) J[ » 2 ' 



x z + u 



2 

(49) 



where x = ojb/j. 

In a similar way the equivalent photon spectrum of a point-like magnetic dipole, 
moving with a constant velocity at a given impact parameter, is calculated in ]56[ |. 
This purely classical result can also be interpreted quantum-mechanically by using 
the corresponding electromagnetic matrix-elements. 

An interesting question is the incoherent contributions due to the protons in- 
side the nucleus. We generalize Eqs. (^7|)-(p9|) in the following way: The static 
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charge distribution is replaced by a transition charge density pfoif) and the Fourier 
transform of it is given by 



f f o(k) = / d d rexpHfcf)p /0 (r), (50) 



where / is some final state of the nucleus and 

z 



ffo(k) = J d 3 r d^}(Oj2S(r-n)^ (OeM-^ (51) 

Z 

= E / d£*/(OexpH*tfi)*o(0, (52) 
i=i 

where £ = ri, T2, • • • , Tz is the set of all proton coordinates (for our argument we 
can neglect the neutrons). We get the total incoherent contribution by summing 
over all states / excluding the ground state f — 0. The sum can be performed using 
the closure relation and we obtain: 

S(k,k') = X)//o(*)//o(*0 



dt\ de £ *5(0 exp(i&*)*/(0*/(0 expi-ik'fyM?) 



h j 



-/oo(fc)/oo(fc') 

= E / * l^o^)! 2 exp^fcrl - ifc'r-) - /* (fc)/ 00 (fc')- (53) 

We split the sum over ij now into one for i = j and one for i ^ j, following pq . 
In the limit of no correlation (not even Pauli correlations) we obtain 

S(k, %') = ZF(k - k') - ZF(k)F{k'), (54) 

where we have introduced a normalized form factor F(k) of the nucleus as F(k) = 

hS)/Z, i.e., F(0) = 1. 

The equivalent photon number due to the incoherent contribution can now be 
written as 

^( W ,6)=«(£) a ^ j M k>JJ' ± S(k, k'), (55) 

where k = (k±,u>/v), k' = (k\ ,ui/v). Using Eqs. ( |47j ) and (p4|) and defining a 
thickness function T z (b) — see Eq. (|13| ) — where the nucleon density ua is now 
replaced by the charge density p), we have 

N lncoh {uj, b) = [ d 2 r 1 _T z {r 1 _)N pomt (uj, b + r ± ) - ZN form (io, b), (56) 



where N^ orm is the usual equivalent photon number for a given form factor (divided 
by Z 2 ), see Eq. (p9|), N pomt denotes the equivalent photon spectrum due to a point 
charge. For 6 = this expression would diverge, and a suitable cut-off has to be in- 
troduced (b m in = Rproton, or another value of bmin-, for which the equivalent photon 
approximation ceases to be valid). If one is integrating over all impact parameters 
(therefore also including central collisions), the total incoherent equivalent photon 
number can be defined: 



incoh 



(«) = Z [n pomt (u;) - n form (uj)] , (57) 
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where n l (ui) — J d 2 bN l (uj,b) and we have used J d 2 r±T z (f±) = Z. 

Incoherent scattering is a well known general phenomenon. E.g., the blue sky is 
caused by the incoherent scattering of sunlight by gas molecules or other randomly 
distributed dipole scatterer (see, e.g., p.418ff of pofl). It can also be formulated 
in a way that makes a connection with the parton model. As an example let us 
look at the contribution from incoherent photon emissions of the quarks in the 
proton. In the parton model a proton in the infinite momentum frame consists of 
partons (quarks, gluons,...). We use the plane wave approach (see where 
the impact parameter dependence is not manifest. The equivalent photon spectrum 
corresponding to a given proton final state can be expressed in terms of the structure 
functions W\ and W%. We only consider the dominant term corresponding to W 2 . 
We are also interested only in inclusive reactions, therefore we integrate over the 
invariant mass of the final state M 2 . One has po| , |58| 

with Q 2 = Q 2 min (ijj) + q\. Since p = p' + q wc obtain 

M 2 = m 2 + 2mv - Q 2 , (59) 

with v = —pq/m. (Since we assume that the photon is emitted, not absorbed, from 
the nuclear system, our sign of q and v is somewhat unconventional.) Introducing 
the scaling variable x = Q 2 /(2mv), we can write in the scaling limit, see e.g. |58fl , 

W 2 {M 2 ,Q 2 ) = -F 2 (x). (60) 
v 



Changing the integration variable dM to dx, we obtain 



a 



1 f^2\Q±\ 



2 



n H = -J d xF 2 (x)-J d Q 2 ^. (61) 

We express the structure function F 2 (x) in terms of the quark distribution func- 
tions f qi \ p (x): 

F 2 (x)=xY / e 2 i f qi \ P (x), (62) 
n 

where is the charge of the quark % (i = u,d, s, • • •). Eq. ( |6l|) now has an intuitive 
interpretation: The proton consist of partons (=quarks) with an momentum fraction 
x, they radiate as pointlike particles. This is described by the dQ 2 integration. 
Taking also the energy loss of the quarks into account in a more complete formula, 
the photon spectrum from a quark is given by 

a 1 + (l-z) 2 , (Q\ 



/W*) = ^ ^^Mt^)' (63) 



'. 2 



where z is the ratio of photon-energy u> and energy of the quark xE. 

According to [ p9[ we choose Q\ to be the maximum value of the momentum 
transfer given by XiX 2 z\z 2 s/4: — m 2 and the choice of the minimum Q\ = 1 GeV 2 
is made such that the photons are sufficiently off shell for the quark-parton model 
to be applicable. The inelastic contribution to the 77-cross section is then given by 

<p 1 = E e ? e J [ 4m2 dx i L 2 dx2 [ im2 dz ^ j im2 dz2 

ij 3 J"l "1"2 ™l , 2 z l 

fq l \p{xi,Q 2 )fq :j \p{xi,Q 2 )fy/ qi (zi)f J /g :j (Z2) 

<7 77 (W 77 = y xix 2 ziz 2 s/4:) (64) 
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(Eq. (3) of |p9[) and where s is the invariant mass of the pp system. Of course 
the photons are now somewhat more off-shell than in the elastic case, and for some 
cases it could be less safe to use the assumption of real photons (q 2 = 0), when 
calculating the cross section for the 77-subprocess a. 




Figure 13: With Q 2 < l/R 2 the photon is emitted coherently from all "partons" 
inside the ion. For Q 2 3> l/R 2 the "partonic" structure of the ion is resolved. 

A similar formula can also be written for the semielastic ( = elastic-inelastic) 
cross section, see ||?||(|. 77-luminosities are calculated according to this procedure 
in both of these papers. The MRSD' parameterization (6l],|6^] for the partonic 
densities is used in ]59[| , in J60[ ] the simple parameterization F2 — 0.161n(l/x) was 
used. One finds that typically the inel-inel contribution is largest, as the charges of 
the partons (quarks) are comparable to the charges of the proton. 

For the proton contribution to the photon spectrum of a heavy nucleus, the 
situation is different. The coherent contribution is proportional to Z 2 , whereas 
they are only proportional to Z for the incoherent one. As Z 1 one expects the 
coherent part to be dominant. 



4 7-A interactions 

The cross section for the collisions of the equivalent photons of one nucleus with 
the other is given by (see Eq. (p5|)): 

f duj 

a = / — n(wW 7 (w). (65) 
J u 

where the equivalent photon number n(u>) is given in Eq. ( pi|) and c 7 (aj) is the 
photonuclear cross section. This gives rise to many interesting phenomena ranging 
from the excitation of discrete nuclear states, giant multipole resonances (especially 
the giant dipole resonance), quasideuteron absorption, nucleon resonance excita- 
tion to the nucleon continuum (see, e.g.. |^,|6^]). Photo-induced processes lead in 
general to a change of the charge-to-mass ratio of the nuclei, and with their large 
cross section they are therefore a serious source of beam loss. Especially the cross 
section for the excitation of the the giant dipole resonance, a collective mode of the 
nucleus, is rather large for the heavy systems (of the order of 100b). For a recent 
discussion see |64). The cross section scales approximately with Z 10 / 3 . (Another 
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serious source of beam loss, the e + e _ bound-free pair creation will be discussed in 
Sec. ^). The contribution of the nucleon resonances (especially the A resonance) 
has also been confirmed experimentally in fixed target experiments with 60 and 200 
GeV/A (heavy ions at CERN, "electromagnetic spallation") J65|-|68|. For details of 
these aspects, we refer the reader to jl6], ^?], [l8[ ^3i| , where scaling laws, as well as 
detailed calculations for individual cases are given. 

Recently the total dissociation cross section for different ion species was studied 
in an experiment at CERN j70). There it was found that this cross section is 
dominated at medium and large Z by the electromagnetic dissociation, with the 
region of the GDR contributing with about 80%. The theoretical calculations, which 
can be considered to be fairly reliable and detailed, see e.g. [fri[ , agree quite well 
with the experiments apart from an additional effect, which can be parameterized 
a s <J a dd = 0.Y1Z barn, a very sizeable effect (For large Z it is even larger than 
the nuclear cross section). It is tempting to guess that the Z dependence is due 
to an incoherent effect of the Z protons in a nucleus. However the corresponding 
incoherent photon number (Eq. (|5^)) is very small for the relevant region b > 
Ri+R.2- Therefore we exclude an incoherent effect as the explanation of the anomaly 
observed in |f70| . 

The interaction of quasireal photons with protons has been studied extensively 
at the electron-proton collider HERA (DESY, Hamburg), with ^/s = 300 GeV 
(E e — 27.5 GeV and E p = 820 GeV in the laboratory system). This is made 
possible by the large flux of quasi-real photons from the electron (positron) beam 
(for a review see pq]). The obtained 7p center-of-mass energies (up to W 1P ~ 
200 GeV) are an order of magnitude larger than those reached by fixed target 
experiments. Similar and more detailed studies will be possible at the relativistic 
heavy ion colliders RHIC and LHC, due to the larger flux of quasireal photons 
from one of the colliding nuclei. In the photon-nucleon subsystem, one can reach 
invariant masses W^n up to WjN,max = V^WmaxEjy « GeV. In the case 

of RHIC ( 197 Au, 7 = 106) this is about 30 GeV, for LHC ( 208 Pb, 7 = 2950) one 
obtains 950 GeV. Thus one can study physics quite similar to the one at HERA, 
with nuclei instead of protons. Photon-nucleon physics includes many aspects, 
like the energy dependence of total cross-sections, diffractive and non-diffractive 
processes (see, e.g., |38f| ). An important subject is elastic vector meson production 
7P — > Vp (with V = p, uj, <j), J/^, ■ ■ •)■ A review of exclusive neutral vector meson 
production is given in [|72| . The diffractive production of vector mesons allows one 
to get insight into the interface between perturbative QCD and hadronic physics. 
Elastic processes (i.e., the proton remains in the ground state) have to be described 
within nonperturbative (and therefore phenomenological) models. It was shown 
in 1 73 1 that diffractive ("elastic") J/^> photoproduction is a probe of the gluon 

density at x « -Mr- (for quasireal photons). Inelastic J/\l/ photoproduction was 

also studied recently at HERA Q. Going to the hard exclusive photoproduction 
of heavy mesons on the other hand, perturbative QCD is applicable. Recent data 
from HERA on the photoproduction of J/^ mesons have shown a rapid increase of 
the total cross section with W^n, as predicted by perturbative QCD. Such studies 
could be extended to photon- nucleus interactions at RHIC, thus complementing the 
HERA studies. Equivalent photon flux factors are large for the heavy ions due to 
coherence. On the other hand, the A-A luminosities are quite low, as compared 
to HERA. Of special interest is the coupling of the photon of one nucleus to the 
Pomeron-field of the other nucleus. Such studies are envisaged for RHIC, sec jis] [Fl| 
where also experimental feasibility studies were performed. 

It is useful to have estimates of the order of magnitude of vector meson produc- 
tion in photon-nucleon processes at RHIC and LHC. Let us assume a cross-section 
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that rises with the jp center of mass energy approximately with a power law: 



f 



(w lN V 

V Wo J 



(66) 



with W chosen to be 1 GeV and « 0.22 for V = p,w, (j), and /3 w 0.8 for V = J/*. 
Also the total hadronic interaction cross section can be parameterized in this form 
with p w 0.16. From Fig. 17 of Q or Fig. 5 of [£EJ one has for cr ss 50 pb for the 
total hadronic cross section, 5 [ib for V — p, 0.5 /ib for V = w, and 10~ 3 ^b for 
y = J/^. Making use of the photon number of Eq. ( p6| ) the total cross section for 
vector meson production on the reaction Z + p — > Z + p+V due to the equivalent 
photon spectrum of the nucleus Z is obtained as 
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(67) 



where some minimum value of the energy of the equivalent photon is used, say 
Wmin = 1 GeV for the total hadronic cross section, u) m in = 2 GeV for V — p, uj 
and (j) an d ^min — 10 GeV for V = J/*, and u) max = Jp/R- For (3 — one obtains 
in a similar way 
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(68) 



The Lorentz factor 7 p of the nucleus Z, as viewed from the proton, is given by 

7p = 2 7 2 -l. (69) 

We assume that a proton and a nucleus Z collide, with the same value 7 (see Fig. |l|, 
where one nucleus is replaced by a proton). We obtain the following numbers, shown 
in Table § 





RHIC 


LHC 


Ion 


I 


Au 


Ca 


Pb 


7 


111 


106 


3750 


2950 


R(fm) 


6 


7 


4 


7 


otot (mb) 


20 


40 


15 


200 


o p (mb) 


2 


3.5 


1.5 


25 


er w , (70 (mb) 


0.2 


0.35 


0.15 


2.5 


cj/* (M b ) 


0.5 


1.5 


3 


40 



Table 2: The expected cross sections for the elastic vector meson production on a 
proton induced by the equivalent photons of a nucleus Z are given for RHIC and 
LHC conditions. Also shown are the total electromagnetic cross sections atot- See 
text for details. 

The numbers in Table |2| refer to the photoproduction on one proton. In AA 
collisions there is incoherent photoproduction on the individual A nuclcons. Shad- 
owing effects will occur in the nuclear environment and it will be interesting to 
study these. There is also the coherent contribution where the nucleus remains in 
the ground state. Due to the large momentum transfer, the total (angle integrated) 
coherent scattering shows an undramatic A 4 / 3 dependence. This is in contrast to, 
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e.g., low energy vk elastic scattering, where the coherence effect leads to an A 2 
dependence, which is relevant for the stellar collapse, see, e.g., | fr6[ , where also a 
pedagogical general discussion of coherence effects is given. In addition there are 
inelastic contributions, where the proton (nucleon) is transformed into some final 
state X during the interaction (see j74|). The experimental possibilities of this at 
RHIC are investigated in 

At the LHC one can extend these processes to much higher invariant masses W, 
therefore much smaller values of x will be probed. Whereas the J/\l/ production at 
HERA was measured up to invariant masses of W ~ 160 GeV, the energies at the 
LHC allow for studies up to « 1 TeV. 

At the FELIX detector at LHC [^| hard diffractive vector meson photoproduc- 
tion can be investigated especially well in AA collisions. In comparison to previous 
experiments, the very large photon luminosity should allow observation of processes 
with quite small jp cross sections, such as T-production. For more details see p3|| . 

C = — 1 vector mesons can be produced in principle by the fusion of three (or, 
less important, five, seven, . . . ) equivalent photons (see Fig. U). The cross section 
scales with Z 6 . This graph has been calculated by I. Ginzburg et al. jnj using the 
methods of [[78|-^0|. It is smaller than the contributions discussed above, even for 
nuclei with large Z. Quite similarly the QED analogue of the C = — 1 mesons, the 
ortho-states of positronium, muonium, or tauonium can be calculated [pl[ . 




Figure 14: By using three (and more) photon processes also states that are forbidden 
in the two-photon fusion process can be produced. Interesting final states f are, for 
example, orthopositronium, orthomuonium, or the vector mesons P,lo,4>. 

One might have thought that these three- and more photon contributions are 
of the same order of magnitude, since Za < 1. But there is another scale: As 
one can see from Eq. ( |76| ) below, there is another factor l/(mb). For electrons 
the important range b w 1/m and this factor therefore is of the order of unity. 
From this one expects that orthopositronium production can be similar in size than 
parapositronium production. For heavier systems the impact parameter range b 
cannot be smaller than R, the nuclear radius. Putting 1/R = A sa 30 MeV (for R « 
7/m), the scale factor is A/M with M the mass of the produced particle M = m M , 
.... This factor is always (much) smaller than one, thus the production via these 
higher order processes is small. Of course the above arguments are rather qualitative 
and should be complemented by more detailed calculations. Such calculations can 
be done using the technique developed for vector meson production in proton-proton 
collisions via gluon-exchange processes [(78|-p0[ . 

As another possibility we mention photon-gluon fusion leading to the production 
of cc and bb quark pairs. It was suggested in Ref. ]82| as a possibility to deduce the 
in-medium gluon distribution. Further studies were done in Refs. ]83|,|84|], and this 
possibility is reviewed in jlfj]. 
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5 Photon-Photon Physics at various invariant 
mass scales 



Scattering of light on light, while absent in classical Maxwell electrodynamics, takes 
place due to quantum effects, like pair creation. At low energies, photon-photon 
scattering is dominated by electron intermediate states, the scattering of light on 
light occurs in higher orders via an electron loop, see, e.g., The lowest order 
process is e + e~ pair creation and is well described by QED. 

Up to now photon-photon scattering has been mainly studied at e + e~ colliders. 
Many reviews |^,^6[|8^] as well as conference reports [|S8|-[9lf exist. The traditional 
range of invariant masses has been the region of mesons, ranging from ir° (m^a = 
135 MeV) up to about r] c (m Vc = 2980 McV). Recently the total 77 — ► hadron cross- 
section has been studied at LEP2 up to an invariant mass range of about 70 GeV 
|92]| . We are concerned here mainly with the invariant mass region relevant for 
RHIC and LHC (see the 77- luminosity figures below). Apart from the production 
of e + e~ (and /i + /i~) pairs, the photons can always be considered as quasireal. The 
cross section section for virtual photons deviates from the one for real photons 
only for Q 2 , which are much larger then the coherence limit Q 2 < 1/R 2 (see also 
the discussion in For real photons general symmetry requirements restrict 

the possible final states, as is well known from the Landau- Yang theorem p^j . 
Especially it is impossible to produce spin 1 final states. In e + e~ annihilation only 
states with J PC = 1~~ can be produced directly. Two photon collisions give access 
to most of the C = +1 mesons. 

The cross section for 77-production in a heavy ion collision factorizes into a 77- 
luminosity function and a cross-section cr 77 (W 77 ) for the reaction of the (quasi)real 
photons 77 — > /, where / is any final state of interest (see Eq. (|27]). When the 
final state is a narrow resonance, the cross-section for its production in two-photon 
collisions is given by 

a yj ^ R (M 2 ) = 8n 2 (2J R + l)r 77 (i?)<5(Af 2 - M 2 R )/M R , (70) 

where J R , M r and T 77 (i?) are the spin, mass and two-photon width of the resonance 
R. This makes it easy to calculate the production cross-section uaa^aa+r of a 
particle in terms of its basic properties. In Fig. [l5| the function ATr 2 dL 17 /dM/M 2 
is plotted for various systems. It can be directly used to calculate the cross-section 
for the production of a resonance R with the formula 



cr AA ^ AA+R [M) = (2J R + ljl 77 • ( 71 ) 



We will now give a general discussion of possible photon-photon physics at rela- 
tivistic heavy ion colliders. Invariant masses up to several GeV can be reached at 
RHIC and up to about 100 GeV at LHC. 

We can divide our discussion into the following two main subsections: Basic QCD 
phenomena in 77-collisions (covering the range of meson, meson-pair production, 
etc.) and 77-collisions as a tool for new physics, especially at very high invariant 
masses. 



5.1 Basic QCD phenomena in 77-collisions 

5.1.1 Hadron spectroscopy: Light quark spectroscopy 

One may say that photon-photon collisions provide an independent view of the 
meson and baryon spectroscopy. They provide powerful information on both the 
flavor and spin/angular momentum internal structure of the mesons. Much has 
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Figure 15: The universal function Ai: 2 dL 11 /dM 11 /M 2 1 is plotted for different ion 
species at LHC (a) and RHIC (b). For the parameters used see table^. 



already been done at e + e~ colliders. For a review see, e.g., |94|. Light quark 
spectroscopy is very well possible at RHIC, benefiting from the high 77-luminosities. 
Detailed feasibility studies exist [fi8|-fjl| . In this study, 77 signals and backgrounds 
from grazing nuclear and beam gas collisions were simulated with both the FRITIOF 
and VENUS Monte Carlo codes. The narrow pj_-spectra of the 77-signals provide a 
good discrimination against the background. The possibilities of the LHC are given 
in the FELIX Lol ||. 

The absence of meson production via 77-fusion is also of great interest for glue- 
ball search. The two-photon width of a resonance is a probe of the charge of its 
constituents, so the magnitude of the two-photon coupling can serve to distinguish 
quark dominated resonances from glue-dominated resonances ("glueballs"). In 77- 
collisions, a glueball can only be produced via the annihilation of a qq pair into a 
pair of gluons, whereas a normal gg-meson can be produced directly, so we estimate 

^(77 ^ M) = T(M -> 77) _ J_ (72) 
(7(77 -> G) T(G -> 77) a 2 ' 

where a s is the strong interaction coupling constant and where G is a "glueball" , M 
a normal qq- meson. Glueballs are produced most easily in gluon-rich environment. 
This is, e.g., the case in radiative J/^f decays, J/ty — ► ^gg. 

In order to form a meson out of the gluon pair, they must first annihilate into a 
qq pair. So we estimate 

r(j/$ -> jG) 1 



r( J/tf -» 7M) a 2 ' 
The "stickiness" of a mesonic state is defined as (see, e.g., 



(73) 



Sx r( j/* - 77) ' (?4) 

We expect the stickiness of all mesons to be comparable, while for glueballs it should 
be enhanced by a factor of about Sg/Sm ~ l/c*s ~ 20, 

In a recent reference |9(| results of the search for /j(2220) production in two- 
photon interactions were presented. There a very small upper limit for the product 
of TjjBk s k s was given, where Bk 3 k 3 denotes the branching fraction of its decay 
into K S K S . From this it was concluded that this is a strong evidence that the 
/j(2220) is a glueball. 
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5.1.2 Heavy Quark Spectroscopy 



For charmonium production, the two-photon width T 77 of r\ c (2960 MeV, J PC — 
is known from experiment. But the two-photon widths of P-wave charmonium 
states have been measured with only modest accuracy. For RHIC the study of r\ c 
is a real challenge JL9|| ; the luminosities are falling and the branching ratios to 
experimental interesting channels are small. 

In Table || (adapted from table 2.6 of ^3|) the two-photon production cross- 
sections for cc and bb mesons in the rapidity range \Y\ < 7 are given. Also given are 
the number of events in a 10 6 sec run with the ion luminosities of 4x 10 30 cm _2 s _1 for 
Ca-Ca and 10 26 cm _2 s _1 for Pb-Pb. Millions of C-even charmonium states will be 
produced in coherent two-photon processes during a standard 10 6 sec heavy ion run 
at the LHC. The detection efficiency of charmonium events has been estimated as 
5% for the forward-backward FELIX geometry |23|], i.e., one can expect detection of 
about 5 x 10 3 charmonium events in Pb-Pb and about 10 6 events in Ca-Ca collisions. 
This is two to three orders of magnitude higher than what is expected during five 
years of LEP200 operation. Further details, also on experimental cuts, backgrounds 
and the possibilities for the study of C-even bottonium states are given in p3]. 



State 


Mass, 


r 77 


a(AA- 


AA + X) 


Events for 10 B sec 




MeV 


keV 


Pb-Pb 


Ca-Ca 


Pb-Pb 


Ca-Ca 


rf 


958 


4.2 


22 mb 


125 fib 


2.2 x 10 7 


5.0 x 10 8 


Vc 


2981 


7.5 


590 nb 


3.8 fib 


5.9 x 10 5 


1.5 x 10 7 


XOc 


3415 


3.3 


160 fib 


1.0 fib 


1.6 x 10 5 


4.0 x 10 6 


X2c 


3556 


0.8 


160 fib 


1.0 fib 


1.6 x 10 5 


4.0 x 10 6 


Vb 


9366 


0.43 


370 nb 


3.0 nb 


370 


12000 


Vob 


9860 


2.5 x 10~ 2 


18 nb 


0.14 nb 


18 


640 


V2b 


9913 


6.7 x 10" 3 


23 nb 


0.19 nb 


23 


76 



Table 3: Production cross sections and event numbers for heavy quarkonia produced 
in a 10 6 sec run in Pb-Pb and Ca-Ca collisions at the LHC with luminosities 1CP 7 
and ^xlO 30 cm~ 2 sec -1 . Adapted from JgJ[/ . 



5.1.3 Vector-meson pair production. Total hadronic cross-section 

There are various mechanisms to produce hadrons in photon-photon collisions. Pho- 
tons can interact as point particles which produce quark- antiquark pairs (jets) (see 
Fig. |l6|a), which subsequently hadronize. Often a quantum fluctuation transforms 
the photon into a vector meson (p,oj,<j), . . . ) (VMD component) opening up all the 
possibilities of hadronic interactions (Fig. |l6|b). In hard scattering, the structure 
of the photon can be resolved into quarks and gluons. Leaving a spectator jet, the 
quarks and gluon contained in the photon will take part in the interaction, some 
examples are given in Figs. |l6| c and d. It is of great interest to study the relative 
amounts of these components and their properties. 

The L3 collaboration recently made a measurement of the total hadron cross- 
section for photon-photon collisions in the interval 5GeV < W 77 < 75GeV J92j. 
It was found that the 77 — >hadrons cross-section is consistent with the universal 
Regge behavior of total hadronic cross-sections. They show a steep decrease in the 
region of low center of mass energy followed by a slow rise at high energies. It is 
parameterized as follows (see Eq. (5) of [p2[ ) 

a tot = A(s/s Q y + B(s/s )-v, (75) 
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(a) 





(c) 



OX. 
g ^V^q 



(d) 



Figure 16: Diagrams showing the contribution to the 77 -^hadron reaction: direct 
mechanism (a), vector meson dominance (b), single (c) and double (d) resolved 
photons. 



with e = 0.0790 ± 0.0011, r) = 0.4678 ± 0.0059, A = 173 ± 7 nb, B = 519 ± 125 nb, 
and so = 1 GeV 2 . Using this together with the effective luminosities (see Table 0), 
one expects about 3 x 10 6 events/GeV at an invariant mass of 25 GeV and still 
4 x 10 4 events/GeV at W 77 = 100 GeV. 

The production of vector meson pairs can well be studied at RHIC with high 
statistics in the GeV region For the possibilities at LHC, we refer the reader 
to 23 and p4j, where also experimental details and simulations are described. 



5.2 77-collisions as a tool for new physics 

The high flux of photons at relativistic heavy ion colliders offers possibilities for the 
search of new physics. This includes the discovery of the Higgs-boson in the 77- 
production channel or new physics beyond the standard model, like supersymmetry 
or compositeness. 

Let us mention here the plans to build an e + e~ linear collider. Such future linear 
colliders will be used for e + e~, ej and 77-collisions (PLC, photon linear collider). 
The photons will be obtained by scattering of laser photons (of eV energy) on 
high energy electrons (« TeV region) (see 97|). Such photons in the TeV energy 
range will be monochromatic and polarized. The physics program at such future 
machines is discussed in |98|], it includes Higgs boson and gauge boson physics and 
the discovery of new particles. 

While the 77 invariant masses which will be reached at RHIC will mainly be 
useful to explore QCD at lower energies, the 77 invariant mass range at LHC - 
up to about 100 GeV — will open up new possibilities. 

A number of calculations have been made for a medium heavy standard model 



Higgs [g9-102|. For masses m# < 2m\y± the Higgs bosons decays dominantly into 
66, whereas a heavier Higgs decays into a W + W~ pair. For the 77 — > H cross 
section we can use Eq. ( |7l|) , where the two-photon width of the Higgs bosons in the 
standard model can be found, e.g., in fl99fl . The calculations, using the integrated 
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luminosity of Table [l], show that for Ca-Ca collisions only about one Higgs boson is 
produced during one year of the LHC operation. Therefore chances of finding the 
standard model Higgs in this case are marginal pi| . 

An alternative scenario with a light Higgs boson was, e.g., given in [103] in the 
framework of the "general two Higgs doublet model". Such a model allows for a 
very light particle in the few GeV region. With a mass of 10 GeV, the 77-width 
is about 0.1 keV (see Fig. 1 of |103| ]). We get 2 x 10 3 events for Ca-Ca collisions, 
40 for p p and 8 for Pb-Pb, with the integrated luminosities of table [l]. The authors 
of |l03| proposed to look for such a light neutral Higgs boson at the proposed low 
energy 77-collider. We want to point out that the LHC Ca-Ca heavy ion mode 
would also be very suitable for such a search. 

One can also speculate about new particles with strong coupling to the 77- 
channel. Large T 77 -widths will directly lead to large 77 production cross-sections, 
see Eq. (|7l|). We quote the papers [104, 105 1. Since the 77-width of a resonance is 
mainly proportional to the wave function at the origin, huge values can be obtained 
for very tightly bound systems. Composite scalar b osons at W 77 ~ 50 GeV are 
expected to have 77- widths of several MeV |104 . 105 1 . The search for such kind of 
resonances in the 77-production channel will be possible at LHC. Production cross- 
section can be directly read off from Eq. (|7|) and Fig. |TJ|. E.g., take W 11 = 50 GeV 
and assume a width of T 77 = 1 MeV, one obtains for a scalar particle ( Jr = 0) 
ccaCa ~ lMeVIO pb keV -1 = 10 nb. With an integrated luminosity of 4 pb -1 in 
the Ca-Ca mode, one obtains 4 x 10 4 events. 

In Refs. [|9[Q 77-processes at pp colliders (LHC) are studied. It is observed 
there that non-strongly interacting supersymmetric particles (sleptons, charginos, 
neutralinos, and charged Higgs bosons) are difficult to detect at the LHC. The 
Drell-Yan and gg-fusion mechanisms yield low production rates for such particles. 
Therefore the possibility of producing such particles in 77 interactions at hadron 
colliders is examined. Since photons can be emitted from protons which do not break 
up in the radiation process (see also Sec. [|) clean events can be generated which 
should compensate for the small number. Formula and graphs for the production 
of supersymmetric particles are also given in [fl6| , where also further references can 
be found. 

In |5^| it was pointed out that at the high luminosity of L = 10 34 cm _2 s _1 at the 
LHC(pp), one expects about 16 minimum bias events per bunch crossing. Even the 
elastic 77 events will therefore not be free of hadronic debris. Clean elastic events 
will be detectable at luminosities below 10 33 cm~ 2 s _1 . This danger of "overlapping 
events" has also to be checked for the heavy ion runs, but it will be much reduced 
due to the lower luminosities. 

Detailed calculations have also been made for the production of a charged 
chargino pair via 77 — * X1X1 ■ The production of these charginos can be stud- 
ied via their decay into a neutralino and a fermion-antifermion pair, Xi ~> Xififj- 
The most clean channel is into muons or electrons. Such an event would therefore be 
characterized by two fermions of opposite charge (e + e _ , p + pT or e ± /i =F ) together 
with an unbalanced transverse momentum. In order to be able to detect the miss- 
ing momentum, a closed geometry is needed. Studies were made for this process 
as a function of the mass of the chargino. In this case the main background — the 
production of a W + W~ pair also decaying into two leptons of opposite charge — 
was studied also. The cross section for this process was found to be 3.6 pb, com- 
parable to the chargino pair production. But the harder momentum distribution of 
this background process can be used to distinguish it from the chargino production. 
Similar calculations have also been made for pp collisions p0[ . 
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6 Diffractive processes as background 



Diffractive processes are an important class of background to 77 final states. As 
the nuclei can remain intact in these collisions, they have the same signature as 
the photon-photon events. Therefore they cannot be distinguished from each other. 
Diffractive events have been studied extensively at HERA for photon-proton colli- 
sions. A future program will also study diffractive processes involving nuclei [106]. 
Diffraction processes in pp and pp are also well known from studies at the Tevatron 
and ISR (CERN). 

Diffractive events at high energies are best described within Regge theory and 
in the language of the Pomeron (see, e.g., 107 1). It is needless to say that the 
possibility to study Photon-Pomeron and also Pomeron-Pomeron collisions in pe- 
ripheral collisions are interesting fields in themselves. Especially photon-Pomeron 
fusion processes could be of interest, as they allow for final states, which are not 
directly possible in photon-photon events (see also the discussion about the diffrac- 
tive vector- meson production in 7 A-collisions in Sec. ^) . Here we restrict ourselves 
to the estimate, how big the contribution of diffractive processes are compared 
to photon-photon events. Of course it is difficult to give quantitative answers at 
present. 

A number of calculations were made within the phenomenological Dual Parton 
Model (DPM) [ 108 1 . These calculations 1 1 OS | have been interpreted, that Pomeron- 
Pomeron fusion dominates over the photon-photon cross section for almost all ions 
used. Only for the very heavy ions, like Pb Pb, the photon-photon process becomes 
comparable. Unfortunately these calculations were made without the constraint 
that the nuclei remain intact in the final state. As the nuclei are only weakly 
bound system and the nuclear interaction strong, it is very likely that a short 
range interaction between them leads to the breakup of the nucleus. More refined 
calculations have been made in the meantime [ 109 1 . The cross sections for diffractive 
processes are then reduced roughly by two orders of magnitude for Ca-Ca and by 
three orders of magnitude for Pb-Pb (both at LHC conditions). Only for proton- 
proton collisions diffractive processes can be expected to dominate over photon- 
photon ones. Particle production from diffractive processes are also studied in 
11C| , 111 I . They also find that the increase of the cross section with mass number 



A is much smaller than that for electromagnetic processes. 

The problem of separating two-photon signals from background has been stud- 
ied in detail for RHIC conditions in [ 1 1 2 1 . Four sources of background have been 
considered: Peripheral (hadronic) nucleus-nucleus collisions, beam-gas interactions, 
gamma-nucleus interactions and cosmic rays. In order to separate signals from back- 
ground, cuts have been developed which utilize the characteristics of two-photon 
interactions; the most important of these cuts arc multiplicity and transverse mo- 
mentum. It was shown in this reference, to which we refer the reader for details, that 
there are high rates of 77-interactions and that the signals can well be separated 
from background. 



7 Electron-Positron Pair production and QED of 
strong fields 

Electrons (positrons) and to some extent also muons have a special status, which 
is due to their small mass. They are therefore produced more easily than other 
heavier particles and in the case of e + e~ pair production lead to new phenomena, 
like multiple pair production. In addition the Compton wave length of the electron 
(~ 386 fm) is much larger than the size of the nuclei (> 7 fm). This also means 
that the virtuality Q 2 of the photons, which ranges from up to the order of l/R 2 
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can be much larger than the electron mass m 2 .. Whereas in all cases discussed 
above we could treat the photons as being quasireal and relate their cross section 
to the photon cross section, this is no longer the case here in general and therefore 
corrections to the EPA are needed. 

The muon has a Compton wavelength of about 2 fm. This length is of the same 
order as the nuclear radius. We therefore expect that the EPA will give more reliable 
results. Both electrons and muons can be produced not only as free particles but 
also into an atomic states bound to one of the ions. 



7.1 Free pair production, Strong field effects and multiple 
pair production 

We consider the e + e~ pair production in the collision of two nuclei with charges Z\ 
and Zi and relative velocity v. For v — > the electrons and positrons can adjust 
adiabatically to the motion of the nuclei and, with sufficiently high charge Z\ and 
i?2 one can enter a supercritical regime, where [Z\ + Z2)a > 1. Such a situation has 
been studied extensively at GSI and late r at Argonne; we refer the reader to the 
(vast) literature on this subject, see [113] where further references are given. We 
here study the opposite region, with v c. A very useful and complete reference 
for this field is p4| . 

The special situation of the electron pairs can already be seen from the formula 
for the impact parameter dependent probability in lowest order. Using EPA one 
obtains |l[] 



PW(6) 



14 (*«)*■ 



1 



9tt 



hi 



( 7<5 



b 2 \ 2m e b 



(76) 



where S ss 0.681 and 7 = 2~/ 2 m — 1 the Lorentz factor in the target frame, one can 
see that at RHIC and LHC energies and for impact parameters of the order of the 
Compton wave length b rs l/m e , this probability exceeds one. It was first described 



in [115] how unitarity can be restored by considering the production of multiple 
pairs[J. 



Multiple pair production was later studied by a number of authors [117-120] 
using different approximations. A general feature found by all was the fact that the 
probability is given to a good approximation by a Poisson distribution: 



P(N,b) 



[P {1) (b)] 
M 



N 



■ exp 



-P«(6) 



(77) 



where P^'Jb) is the single pair creation probability from perturbation theory, see, 
e.g., Eq. (|76|). Deviations from this Poisson form were studied in [12C] and were 
found to be small at high energies. 

The impact parameter dependence of the lowest order process was calculated 
in [121, 122] (see also Fig. |l7|), the total cross section for the one-pair production 
in [123], for one and multiple pair production in [124] (see Fig. [j^). Of course 



the total cross section is dominated by the single pair production as the main 
contribution to the cross section comes from very large impact parameters b. On 
the other hand one can see that for impact parameters b of about 2R the number 
of electron-positron pairs produced in each ion collision is about 5 (2) for LHC 
with Z = 82 (RHIC with Z = 79). This means that each photon-photon event — 
especially those at a high invariant mass — which occur predominantly at impact 
parameters close to b > 2R — is accompanied by the production of several (low- 
energy) e + e~ pairs. 



1 It is interesting to remark that the fact that pair production in ion collisions grows beyond 
the unitarity limit was already observed by Heitler (lift. Of course at that time available energies 
made this only "of academic interest" . 
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Figure 17: The impact parameter dependent probability to produce N e + e~ -pairs 
(N = 1,2,3, A) in one collision is shown for both RHIC (a,j = 106, Au-Au) and 
LHC (b,-y = 2950, Pb-Pb). Also shown is the total probability to produce at least one 
e + e~-pair. One sees that at small impact parameters multiple pair production can 
be dominant over single pair production. 
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Figure 18: The total cross sections for the multiple pair production of up to three 
pairs is shown as a function of the Lorentz factor 7. Shown are the results for a 
Pb-Pb collision. 
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Figure 19: The single differential cross section for e + e~ pair production at LHC 
energies. Shown are the cross sections as a function of the energy of either electron 
or positron (a), and as a function of the angle with the beam axis (b). 



As the total cross section for this process is huge (about 200 kb for Pb at LHC, 
30 kb for Au at RHIC), one has to take this process into account as a possible 
background process. Most of the particles are produced at low invariant masses 
(below 10 MeV) and into the very forward direction (see Fig. |l9|). Therefore most 
of them leave the detector along the beam pipe and are not observed. On the other 
hand a substantial amount of them is left also at high energies, e.g., above 1 GeV. 
These QED pairs also constitute a potential hazard for detectors at the colliders. 
In Table |] we show the cross section for this process with the energy of either the 
electron or the positron above a certain threshold. Singles angular distributions 
of electrons (positrons) are calculated for peripheral collisions using EPA in [125]. 
Numerical results in the relevant energy and angular ranges are presented there. 
The physics is discussed in terms of easy to handle analytical formulae. 



E thr (GeV) 


cr(Pb-Pb,LHC) 


cr(Ca-Ca,LHC) 


0.25 


3.5 kb 


12 b 


0.50 


1.5 kb 


5.5 b 


1.0 


0.5 kb 


1.8 b 


2.5 


0.08 kb 


0.3 b 


5.0 


0.03 kb 


0.1 b 



Table 4: The cross sections of e + e~ pair production with both electron and positron 
having an energy above a certain threshold value. Shown are results for both LHC 
for two different ion species. 



Differential production probabilities for 77-dileptons in central relativistic heavy 
ion collisions are calculated using EPA and an impact parameter formulation and 
compared to Drell- Yan and thermal ones in Ji5,126,127|. The very low p± values and 
the angular distribution of the pairs give a handle for their discrimination. Nuclear 
stopping leads to bremsstrahlung pair-production and some modification of the 77 
dilepton spectra. For details we refer the reader to these references. In [ 128 1 the low 
energy dilepton spectrum in 16 GeV 7r-p collisions was studied using the two-photon 
mechanism. It was found that this mechanism could not explain the experimental 
data [12S]. This is in contrast to the findings of 130 . 
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In the Bethe-Heitler process 7 + Z — > e + e~ + Z higher-order effects are well 
known. Using Sommerfeld-Maue wave functions higher order effects are taken into 
account. This leads to a modification of the Born result. E.g., the total cross section 
(no screening) for to ^> m e is given by (see [ 13 1 1 ) 
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with the higher-order term given by 



f(Za) 



(78) 



f(Za) = (Za 



^ n(n 2 + (Za) 2 ) 



(79) 



and r e = a/m e is the classical electron radius. As far as total cross sections are 
concerned the higher-order contributions tend to a constant. 

Using those results a modification of e + e~ pair creation in Zi + Z 2 collisions 
with respect to the lowest order result was obtained Q. Such a treatment was 
not symmetrical with respect to Z\ and Z 2 and an ad hoc symmetrization was 
introduced there (see Eq. (7.3.7) of 0). 

A systematic way to take leading terms of higher order effects into account in 



e + e pair production is pursued in [132] using Sudakov variables and the impact- 
factor representation. 



Nonperturbative effects are studied also in a light-front approach [ 133 1 . In this 
approach a gauge transformation on the time-dependent Dirac equation is per- 
formed, in order to remove the explicit dependence on the long-range part of the 
interaction. Similar approaches are also studied in [134. 135 1. Numerical evaluation 



of the non-perturbative effects will be considered in a future work. 



In this context the paper [136] is of interest. In this work the electromagnetic 
field of a particle with velocity v is calculated, see e.g., the textbook result of p0| . 
Then the limit v — > c is performed. This corresponds to the electromagnetic fields 
of a massless particle, which can be regarded as an "electromagnetic shock wave" . 
The results of this paper are very much reminiscent of the sudden approximation 
in the semiclassical theory; for a connection of semiclassical and eikonal methods 
see also 1 137 . 



7.2 Bound-Free Pair Production 

The bound-free pair production, also known as electron-pair production with cap- 
ture, is a process, which is also of practical importance in the collider. It is the 
process, where a pair is produced but with the electron not as a free particle, but 
into an atomic bound state of one of the nuclei. As this changes the charge state 
of the nucleus, it is lost from the beam. Together with the electromagnetic dissoci- 
ation of the nuclei (see Sec. f|) these two processes are the dominant loss processes 
for heavy ion colliders. 

In [nj an approximate value for this cross section is given as 



10 1 z e cxp(27rZ 2 a) - 1 



In ( 7( 5/2) -jj 



(80) 



where only capture to the if -shell is included. The cross section for all higher shells 
is expected to be of the order of 20% of this cross section (see Eqs 7.6.23 and 24 
of |). 

The cross section in Eq. (|8(]) is of the form 

a = A\nj + B. (81) 
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This form has been found to be a universal one at sufficient high values of 7. The 
constant A and B then only depend on the type of the target. 

The above cross section was found making use of the EPA and also using ap- 
proximate wave function for bound state and continuum. More precise calculations 
exist [138 -143 1 in the literature. Recent calculations within PWBA for high values 
of 7 have shown that the exact first order results do not differ significantly from 
EPA results [144, 145]. Parameterizations for A and B [ 139| , 141 for typical cases 
are given m Table§ 



Ion 


A 


B 


a(Au,~f = 106) 


a(Pb,j = 2950) 


Pb 


15.4b 


-39.0b 


115 b 


222 b 


Au 


12.1b 


-30.7b 


90 b 


173 b 


Ca 


1.95mb 


-5.19mb 


14 mb 


27.8 mb 





4.50^b 


-12.0/ib 


32 nb 


64.3 ^b 



Table 5: Parameters A and B (see Eq. j\8J\)) as well as total cross sections for the 
bound-free pair production for RHIC and LHC. The parameters are taken from \14aI- 



For a long time the effect of higher order and nonperturbative processes has 
been under investigation. At lower energies, in the region of few GeV per nucleon, 
coupled channel calculations have indicated for a long time, that these give large 
contributions, especially at small impact parameters. Newer calculation tend to 
predict considerably smaller values, of the order of the first order result and in a 
recent article Baltz [146] finds in the limit 7—5-00 that contributions from higher 
orders are even slightly smaller than the first order results. 

The bound-free pair production was measured in two recent experiments at the 
SPS, at 7 = 168 (147) and at 7 w 2 Jl4|,|l4|). Both experiments found good 
agreement between measurement and calculations. 

A similar process was recently used at LEAR (CERN) to produce antihydrogen. 
An antiproton beam with a momentum of 1.94 GeV/c hit a Xenon target (Z — 54) to 
produce and detect antihydrogen in the bound free pair production mechanism ]150| . 



The same technique is also used at Fermilab [151], where it is also planned to 
measure the Lamb shift in antihydrogen as a test of CPT invariance [152,153]. 
Their results are in good agreement with the recent calculations [144, 145| 



We note that electron and positron can also form a bound state, positronium. 
This is in analogy to the 77-production of mesons (qq states) discussed in Sec. [|. 
With the known width of the parapositronium r((e + e~) n= i So — > 77) = mc 2 a 5 /2, 
the photon-photon production of this bound state was calculated in [ 1 54 1 . The 
production of orthopositronium, n = l 3 Si was calculated recently [j77f . 

As discussed in Sec. |^ the production of orthopositronium is only suppressed 
by the factor (Za) 2 , which is not very small. Therefore one expects that both 
kind of positronium are produced in similar numbers. Detailed calculation show 
that the three-photon process is indeed not much smaller than the two-photon 
process @||. 



8 Conclusion 

In this review the basic properties of electromagnetic processes in very peripheral 
hadron-hadron collisions (we deal mainly with nucleus-nucleus collisions, but proton 
proton collisions are also treated) are described. The method of equivalent photons 
is a well established tool to describe these kind of reactions. Reliable results of 
quasireal photon fluxes and 77-luminosities are available. Unlike electrons and 
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positrons heavy ions are particles with an internal structure. We have described 
how to treat effects arising from this structure, and we conclude that such effects 
are well under control. A problem, which is difficult to judge quantitatively at 
the moment, is the influence of strong interactions in grazing collisions, i.e., effects 
arising from the nuclear stratosphere and Pomeron interactions. 

The high photon fluxes open up possibilities for photon-photon as well as photon- 
nucleus interaction studies up to energies hitherto unexplored at the forthcoming 
colliders RHIC and LHC. Interesting physics can be explored at the high invariant 
77-masses, where detecting new particles could be within range. Also very interest- 
ing studies within the standard model, i.e., mainly QCD studies will be possible. 
This ranges from the study of the total 77-cross section into hadronic final states 
up to invariant masses of about 100 GeV to the spectroscopy of light and heavy 
mesons. 

We also reviewed dilepton production in very peripheral collisions. This is es- 
sentially well understood and gives rise to an experimental background. Multiple 
pair production is a strong field effect of principle interest. Pair production with 
capture is, in addition to nuclear fragmentation (sometimes called "Weizsackcr- 
Williams process"), a source of beam loss in the collider operation. 

The future is coming soon. RHIC will be operational in only one year, LHC in 
approximately seven years. Therefore the planning of the experiments and necessary 
detectors for this kind of physics has to be done now. With the new data and 
new insights, that will come from these experiments, new work and theoretical 
understanding will be required. As an ancient motto says: "No surprise would be 
a surprise" . 
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